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ABSTRACT 

The outer halo regions of early-type galaxies carry key information about their past 
accretion history. However, spectroscopically probing the stellar component at such 
galactocentric radii is still challenging. Using Keck/DEIMOS, we have been able to 
measure the metallicities of the stellar and globular cluster components in 12 early-type 
galaxies out to more than 10 Rj,. We find similar metallicity gradients for the metal- 
poor and metal-rich globular cluster subpopulations, suggesting a common formation 
process for the two subpopulations. This is in conflict with most current theoretical 
predictions, where the metal-poor globular clusters are thought to be purely accreted 
and metal-rich globular clusters mostly formed in-situ. Moreover, we find that the 
globular cluster metallicity gradients show a trend with galaxy mass, being steeper in 
lower-mass galaxies than in higher-mass galaxies. This is similar to what we find for the 
outermost galaxy stars and suggests a more active accretion history, with a larger role 
played by major mergers, in the most massive galaxies. This conclusion is qualitatively 
consistent with expectations from two-phase galaxy assembly models. Finally, we find 
that the small difference in metallicity between galaxy stars and metal-rich globular 
clusters at 1 Rc may correlate with galaxy mass. The origin of this difference is not 
currently clear. 

Key words: galaxies: abundances - galaxies: elliptical and lenticular, cD - galax¬ 
ies: formation - galaxies: evolution - galaxies: star clusters: general - galaxies: stellar 
content. 


1 INTRODUCTION 

The currently most popular view of massive early-type 
galaxy (ETG) formation is the “two-phase” scenario. The 
first phase consists of an early (2 > 2) dissipative collapse of 
gas that forms the bulk of present-day galaxy stars. The for¬ 
mation of these “in-situ” formed stars is driven by gas-rich 


mergers and/or the gravitational in-fall of gas ( 

De Lucia & 

Blaizotl 20071 |Dekel, Sari & Ceverino |2009 

Zolotov et al. 

2009| IKhochfar & Silk|2009| |Font et a .||2011| |Cooper et al. 
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20151. The second phase, by contrast, involves the accretion 


of stars formed “ex-situ” (i.e. in external satellites), gener¬ 
ally during gas-poor minor merger events. Simulations have 
found that such stars are mostly accreted in the outer re- 


gions of the galaxy at 2 

< 2 (Naab, Johansson & Ostriker 

20091 |Oser et al.|20101|Font ei 

al.|20111 ILackner et al.|2012| 

Navarro-Gonzalez et al. 

2013 

1. In massive galaxies, an ac- 


cretion of stars during major mergers events is predicted 
to strongly reduce the stellar specific angular momentum 


Khochfar et al. (20111. As a consequence, massive slow ro¬ 


tators are thought to have experienced more major mergers 
than the fast rotators with similar mass. 
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2 Pastorello et al. 


We acknowledge that the separation of galaxy formation 
into two phases might be somewhat artificial and too sim¬ 
plistic. This is because mergers happen continuously, with 
the merger rate and gas fraction decreasing steadily in time. 
Therefore, the separation of galaxy formation in two phases 
at z « 2 may be seen as the epoch that roughly denotes a 
shift from frequent gas-rich accretion episodes to more dis¬ 
crete gas-poor merger events. 

While the study of the bright central regions of ETGs 
has provided strong constraints on their early formation 
phase, only the exploration of the outer regions of galax¬ 
ies will help to understand the second accretion phase. In 
addition, more than 90 per cent of the total mass and angu¬ 
lar momentum, and half of the stellar mass lie outside the 
half-light radius (Ro). Therefore, most of the information 
about a galaxy’s formation and evolution is encoded in its 
outer regions. 

The two phase formation scenario makes prediction for 
the stellar metallicity ([.Z/H]) gradients in ETGs. The dif¬ 
ferent phases of formation are expected to produce different 
metallicity profiles for the “in-situ” and “ex-situ” compo¬ 
nents. For instance, the first dissipative phase is supposed 
to create a metal-rich “in-situ” component in the inner re¬ 
gions, where the density is higher and the stellar chemical 


evolution is faster (Chiosi & Garraro 2002 Zolotov et al. 


2010). This “in-situ” component will also reveal a steep ra¬ 


dial metallicity gradient, because of the rapidly decreasing 


local density with galactocentric radius (Kobayashi 2004 
Pipino et al.|2010). 


On the other hand, the “ex-situ” stars form mostly 
in lower-mass, relatively metal-poor, external systems. In 
present-day ETGs, these stars contribute to lowering the 
metallicity in the regions where they have been accreted, 
i.e., at large galactocentric radii. These outermost regions 
are expected to show generally shallow metallicity gradients 
because of the dynamical mixing of pre-existing and accreted 
stars I Kobayashil2004 Di Matteo et al.|2009 Cooper et al. 


2010 Font et al. 20111. 


In galaxies that experienced accretion, a transition be¬ 
tween the “in-situ” metal-rich and the “ex-situ” metal-poor 
dominated regions is expected. The radial position of such 
a transition, which marks where the accreted component 
starts dominate, is still uncertain. For example, the simu¬ 


lations of Hirschmann et al. (2015) show that this transi¬ 
tion radius could lie from 4 to more than 8 Re. Because of 
their different accretion histories, the metallicity profiles of 
low mass galaxies are expected to be systematically different 
from those of their high mass counterparts. 

Finally, major merger events in massive galaxies will 


weaken any pre-existing stellar metallicity gradient (IWhite 

1978 

|1980||Kobayashi||2004||Ui Matteo et al.||2009||Rupke, 

Kew 

ey & Barnesj 2010| Navarro-Gonzalez et al. 2013| 

Hirschmann et al.|2015 

1 . This is because as the mass of the 


two galaxies involved in the merger are similar, their stars 
will also show similar metallicities but will be mixed at all 
galactocentric radii in the final galaxy. 

To date, just a handful of studies have spectroscopi¬ 
cally probed the stellar metallicity at large radii, and most of 
them just in the high-mass galaxy range (i.e. M* > 10^^ M©, 
Goccato, Gerhard fc Arnaboldi|2010[[Greene et aI.|2012||Mi 


hosetakl 20131. An exception to this is Pastorello et al. 


{1 < R ^ 2.5Ro) has been spectroscopically measured also 
in lower-mass galaxies (i.e. with stellar mass down to almost 
10^° Mq). In this work, the outer metallicity gradient was 
found to correlate with galaxy mass, becoming steeper in 
low-mass galaxies. This may confirm that high-mass galax¬ 
ies accrete more, and higher-mass, satellites than their low- 
mass counterparts. However, the study of [Pastorello et al.| 
(2014) did not reach the galactocentric radii where the ac¬ 


creted component is expected to dominate. 

For a small number of very nearby galaxies, stellar 
metallicity studies at very large radii have been based on 
the photometric analysis of resolved red-giant branch (RGB) 
stars (Elson 1997 Harris et al. 2007a|b Bird et al. 2015| 


Rejkuba et al.||2014 Peacock et al.||2015 l. A transition be¬ 

tween possible “in-situ” and “ex-situ” stars has been no¬ 
ticed around 12 Re in NGC 3379 ( Harris et al.|2007b ). In a 
spectroscopic study of M31 RGB stars, Kalirai et al. (2006) 


found a transition in the stellar metallicity at « 30 kpc from 
the galaxy centre. As this method requires resolved stars, it 
is limited to the local Universe, and just a handful of galaxies 
have been studied in this way. 

A way to overcome such issues and explore a higher 
number of galaxies out to large radii is to use a different 
tracer of the galaxy stellar populations. In particular, glob¬ 
ular clusters (GGs) are a viable way to explore the metal¬ 
licity at large galactocentric radii, since they are bright and 
dense stellar systems. For this reason, their study can extend 
beyond the detectable galaxy stellar light, and their metal¬ 


licities can be measured out to more than 10 Re (Usher et al. 


2012). Since GGs are generally very old (i.e. age > 10 Gyr), 


they have survived the evolutionary processes that formed 
their present-day host galaxies. 

Almost all the GC systems of galaxies with M* > 
10^° Mq present a bimodality in colour and metallicity, with 
blue GCs (bGCs) showing lower metallicities than red GGs 
(rGG s) ([Strader et al.|[2006| [Brodie Sz Strader||2006| |Peng| 
|et al.|[^006[ and references therein). The origin of this bi¬ 
modality is still unclear, but it is commonly believed that 
it is intrinsically linked with the host galaxy’s formation 
history (e.g. Strader et al.|2005 Brodie fc Strader|2006 1. In 


general, rGGs seem to be linked with the host galaxy bulge’s 
stellar population, having similar spatial distribution, kine¬ 
matics and metallicity, while bGCs seem to share several 


properties with the galaxy halo component (e.g., Faifer et al. 

2011 

Strader et al.||2011||Arnold et al.||2011||Forbes et al. 

2011 

Forbes, Ponman & O’Sullivan||2012| Pota et al.||2013| 

Peacock et al.|2015 

Kartha et al.||2014 

Pota et al.|2015 1. 


(2014), where the stellar metallicity of the outer region 


In this paper we explore galaxy stellar populations at 
large radii using both stellar and GC metallicities. In partic¬ 
ular, we find that the steepness of the GC metallicity radial 
profile is higher in low mass galaxies, qualitatively similar 
to that between stars and galaxy mass found by jPastorelloj 
|et al.| ( |2014[ ). Furthermore, we measure the metallicity off¬ 
set between the galaxy stars in the innermost regions and 
the rGGs. We find that this offset correlates with the galaxy 
mass, in the sense that higher-mass galaxies present larger 
offsets. 

The paper is structured as follows. In Section we 
present the data and the methods we use to obtain the 
metallicity of both stars and GCs. In particular, in Section 
|2.6| we compare our spectroscopic metallicities with photo¬ 
metric metallicity profiles from the literature in order to 
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Stellar and globular cluster metallicities 3 


assess the possible selection bias for our sample of spectro¬ 
scopic targets. Section [^presents the results we obtain from 
the radial metallicity profiles for both the stellar and the GC 
components. In Sectionj^we present our findings and discuss 
how they compare with the current theoretical predictions 
and literature results. Finally, in Sectionj^we provide a sum¬ 
mary of the results. In addition, in Appendix A we present 
an updated metallicity map for the galaxy NGC 3115, re¬ 
flecting new extra data points obtained since the publication 


of the original in Pastorello et al. (20141. 


2 DATA ACQUISITION AND ANALYSIS 
2.1 The SLUGGS sample 

In this work we present the metallicity radial profiles of both 
the stellar and the GC components for 12 galaxies from the 
ongoing SAGES Le gacy Unifying Glob ulars and GalaxieS 
(SLUGG^ survey ( Brodie et al.]2014 1. The survey sample 
includes 25 early-type galaxies (i.e. ellipticals and lenticu- 
lars) within a distance D < 30 Mpc. One of the aims of this 
survey is the study of the stellar populations in both GCs 
and field stars. We select the sub-set of 12 SLUGGS galax¬ 
ies for which the metallicities of at least 15 spectroscopically 
confirmed GCs have been measured. 

The SLUGGS dataset presented in this paper has been 
obtained over 9 years and 25 observing runs. The GC pho¬ 
tometry is from Usher et al. ( 2012[) f or all the galaxies, 
except NGC 4473 (Alabi et al. 20151, NGC 4649 (Pota 
et al. 20151 and NGC 1023 (this work). Table 0 presents 
the adopted parameters for the galaxies in this paper. 


2.2 DEIMOS spectra 

In order to reduce the DEIMOS spectra, we adopt a modi 
fied version of the IDL spec2D pipeline (Cooper et al.|2(il2 


[Newman et al.||2013[ ), as presented in [Arnold et al.| ( |2014p . 
In multi-slit mode, from a single DEIMOS slit, it is pos¬ 
sible to obtain the spectra of both the target object and 
the background light. While the former is used to study the 
GC stellar population properties, the latter includes light 
from both the galaxy’s stars and the night sky. Adopting 
the “Stellar Kinematics with Multiple Slits” (SKiMS) tech- 


nique (Norris et al. 

2008 

Proctor et al. 

2009 Foster et al. 

2009 

[Pastorello et al.|2014| [Arnold et al.|2014| [Foster et al. 

2015 

I, it is possible to separate the two components in the 


background spectra and obtain an individual spectrum of 
the GC and the galaxy starlight at each slit location. 


ATLAS^® metallicity data for NGC 4365 kindly provided by 
Harald Kuntschner (Kuntschner et ah, in prep.). When the 
galaxy metallicity is not available from the SAURON and 
ATLAS^° surveys, we adopted the 
empirical correction relation: 


Pastorello et al. 


(20141 


A[^/H] = 3.78 • log CTO - 8.08 


( 1 ) 


where the values for the central velocity dispersion ao are 
those in Table [T] 

From these metallicity values, we obtained 2D maps of 


the galaxy metallicity using the kriging technique (see Pas- 
torello et al.[|2014 for details). Kriging maps take into con¬ 
sideration the non-optimal spatial sampling of the available 
datapoints. In this way we are able to obtain pseudo-IFU 
maps out to several effective radii and derive radial metal¬ 
licity profiles. For this purpose, we are adopting the photo¬ 
metric parameters in Table 

Throughout the paper, the galactocentric radius is ex¬ 
pressed in terms of circularized radius, defined as: 


R= \ l qx^ + — 


( 2 ) 


where q is the axial ratio and x, y are the coordinates along 
the major and the minor axes of the galaxy, respectively. 
For consistency, the same definition of radius is also used 
to convert the galactocentric radii of literature metallicity 
values, when these are given on different scales. 


2.4 GC metallicities 

The metallicities of the GCs presented in this work are 
mostly those presented in Usher et al. (20121, Usher et al. 


(20131 and Usher et al. (20151. These are supplemented 


with new data obtained from recent Keck observations re¬ 
duced and analysed using the methods described in [Usher[ 
et al. (20121. In particular, we add new data for NGC 1023, 


NGC 1407 and NGC 3115. From the total sample, we select 
the galaxies with more than 15 spectroscopically confirmed 
GCs with S/N > 10 A in the wavelength range from 8400 
to 8500 A. 

In order to separate the blue and the red GC subpop¬ 
ulations we adopt a split in [g — i) colours as defined in 
Table [T] Such colour values are obtained from the best fit¬ 
ting bimodal colour distributions calculated with the Gaus¬ 


sian Mixture Modelling (GMM) code of Muratov & Gnedin 
( [2010[ ), as in [Usher et al.[ ( |2012| ). The whole sample of GC 
CaT and metallicities used in this paper are presented in 
Table Hi 


2.3 Stellar metallicities 

The SKiMS metallicity data obtained from the CaT index 
have been corrected following the prescriptions in section 


3.4 of Pastorello et al. (20141. In particular, in the galax¬ 
ies for which SAURON or ATLAS®° data are available, we 
match the radial profiles within 1 Re in order to minimize 
the metallicity offset. The SAURON metallicity data are the 


same used in Pastorello et al. (20141, to which we add the 


http://sluggs.swin.edu.au 


2.5 Stellar and GC metallicity scales 

The stellar and the GC metallicities used in this work 
are defined on slightly different metallicity scales. While 
SAURON/ATLAS^® metallicities are obtained from Lick in¬ 
dices and stellar population models from Schiavon (20071, 


the Pastorello et al. (20141 stellar metallicities are from the 


CaT index and the SSP models of Vazdekis et al. (20031. 
In addition, an empirical correction is applied to the [Pas-] 
torello et al. metallicities to correct for an offset with 

SAURON/ATLAS^® metallicities in the overlapping radial 
regions, which is likely driven by the IMF. GC metallicities 
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4 Pastorello et al. 


Galaxy 

(1) 

Re 

(arcsec) 

(2) 

Type 

(3) 

Dist. 

(Mpc) 

(4) 

^kpc 

(km s“^) 

(5) 

Mg 

(mag) 

(6) 

Zog(M*/M0) 

(dex) 

(7) 

(ff-*) 

(mag) 

(8) 

NGC 1023 

48 

SO 

11.1 

183 

-20.15 

10.92 

0.96 

NGC 1400 

28 

El/so 

26.8 

236 

-20.36 

11.03 

0.95 

NGC 1407 

63 

EO 

26.8 

252 

-21.61 

11.50 

0.98 

NGC 2768 

63 

E6/S0 

21.8 

206 

-21.02 

11.20 

0.95 

NGC 3115 

35 

SO 

9.4 

248 

-20.10 

10.91 

0.93 

NGC 3377 

36 

E5-6 

10.9 

135 

-19.02 

10.42 

0.93 

NGC 4278 

32 

El-2 

15.6 

228 

-19.99 

10.83 

0.86 

NGC 4365 

53 

E3 

23.1 

253 

-21.30 

11.40 

0.91 

NGC 4473 

27 

E5 

15.2 

189 

-19.82 

10.82 

0.99 

NGC 4494 

49 

El-2 

16.6 

157 

-20.45 

10.96 

0.99 

NGC 4649 

66 

E2/S0 

16.5 

308 

-21.32 

11.46 

0.99 

NGC 5846 

59 

EO-l/SO 

24.2 

231 

-21.02 

11.31 

0.96 


Table 1. Galaxy parameters. The columns present: (1) Galaxy name. (2) Effective radius. (3) Morphological type, (4) Distance and (5) 
central velocity dispersion from |Brodie et ar] < |201^ . (6) Total B-band magnitude from |de Vaucouleurs et al.| l |199l] |, corrected for galaxy 
extinction. (7) Total logarithmic stellar mass obtained from the X-band extinction corrected absolute magnitude from the 2MASS source 
catalo gue l|Jarrett et aL]|2000[ |, obtained following the same approach as in [Pastorello et al. ||2014[|, assu ming the distances in [Brodiej 
^SilJ (|20^1^||. (8) Colour marking the split between the blue and red GC subpopulations from[TJs ler et al.j l|2012||, except for NGC 1023, 
NGC 4473 and NGC 4649 (calculated in this work). 


are also obtained from the CaT index and the SSP models 
of Vazdekis et al. (20031 (for details, see Usher et al. 

Any differences in the metallicity scale between stellar and 
GC metallicities may, in principle, lead to systematic offsets 
between the two components. 


2.6 Comparison between spectroscopic and 

photometric GC metallicity radial profiles 

Our sample of GC spectroscopic metallicity may suffer from 
a selection bias at low metallicities, given that this selection 
is based on the detection of metallicity-dependent absorp¬ 
tion lines. However, the GC metallicity of several galaxies 
in our sample has already been studied photometrically. For 
these galaxies we compare our spectroscopic GC metallici¬ 
ties with the photometric GC metallicities from the litera¬ 
ture, in order to check whether or not our sample is repre¬ 
sentative of the whole GC population. Most GCs are pre¬ 
dominately old, so photometric colours are a good proxy for 


metallicity (Puzia et al. 2005 Strader et al. 2005 Spitler 


et al. 20081, although the relationship is non-linear and 
seems to vary galaxy-to-galaxy ( Usher et al.|2015 |. 

In the case of NGC 1407, Forbes et al. (20111 pre¬ 


sented {g — i) radial trends of the red and blue GCs from 
Subaru/Suprime-Cam. We convert these colours into metal¬ 
licities following equation 10 in Usher et al. (20121. In their 


tract these profiles from their figure 3 and convert them 
into metallicity profiles using equation 10 in [Usher et al.| 
( |2012[ ). jJennings et al.| ( |2014[ ) presented the (g — z) colour 
profiles obtained from HST /ACS photometry. To convert 
these colours into total metallicities, [^/H], we adopt equa¬ 
tion B1 in Usher et al. (|2012[). 


Usher et al. (20131 presented the colour gradients of 


NGC 4278 rGCs and bGCs from both HST/ACS and 
Subaru/Suprime-Cam photometry. From these, we obtain 
the (g — z) and {V — I) colour gradients for both the GC sub¬ 
populations, which we convert into total metallicities [^/H] 
with equation B1 and B3 in Usher et al. (20121, respectively. 


The case of NGC 4365 is peculiar, since this galaxy 


hosts three, rather than two, GC subpopulations (Puzia 
|et al.|[?002[ ). |Blom, Spitler fc Forbes] ( |2012[ ) presented the 
(g — i) colour profiles obtained from Subaru/Suprime-Cam 
photometry for the three subpopulations. For the purposes 
of a spectroscopic/photometric metallicity comparison, we 
consider only the rGC and the bGC colour profiles in their 
work and exclude the central “green” GCs. We extract these 
colour profiles from their figure 16 and convert them into 


metallicity profiles with equation 10 in Usher et al. (20121. 


Similarly, Pota et al. (20151 present the {g — z) colour 


profiles for both the rGC and bGC subpopulations in 
NGC 4649. These colours are obtained from HST/ACS, 
Subaru/Suprime-Cam and CFHT/MegaCam photometry. 
As per Jennings et al. (20141, we convert Pota et al.l (2015 \ 


work, Forbes et al.j (20111 found that the colour gradients colours into total metallicities \Z/'H\ using equation B1 in 


of both rGCs and bGCs are steeper for i? < 10 arcmin and 
flatten in the outer regions. They were concerned that the 
flattening of the rGC colour gradient in the outer regions 
may be a consequence of strong contamination (more than 
50 per cent) outside of 10 arcmin from the galaxy centre. We 
reduce the probability of this case, since we spectroscopically 
confirm the presence of more than few rGCs out to almost 
15 arcmin in this galaxy. 

The two subpopulations of GCs in NGC 3115 have been 


studied by Arnold et al. (20111 and Jennings et al. (20141. 
Arnold et al. (20111 presented the {g — i) colour profiles 
obtained from Suprime-Cam/Subaru photometry. We ex¬ 


Usher et al. (2012 I. 

In Figure the comparison between photometric and 
spectroscopic metallicities is presented for NGC 1407, 
NGC 3115, NGC 4278, NGC 4365 and NGC 4649. In all 
the galaxies, the bGC average spectroscopic metallicities are 
generally consistent with the photometric metallicities. We 
therefore exclude the possibility of a significant bias at low- 
metallicity for our spectroscopic bGC sample. However, the 
spectroscopic rGC metallicities are generally lower than the 
photometric ones at the level of « 0.3 dex. This may be a 
consequence of (1) non universal colour-metallicity relations 
for GCs, as discussed in Usher et al. (20151, (2) the spectro- 
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Figure 1. GC radial metallicity profiles. Red circles and blue crosses show the spectroscopic metallicities for the rGCs and bGGs, 
respectively. The blue and red points with error bars show the average uncertainties on the rGG and bGC spectroscopic metallicities, 
respectively. All the points are plotted at their elliptical galactocentric radius scaled by Re. In all the panels the red and blue lines show 
the literature photometric metallicity profiles of rGCs and bGCs, respectively. 


scopic sample not being representative of the GC subpop¬ 
ulation, and/or (3) an imperfect classification of GCs into 
the two subpopulations (which is based on the GC colours). 
In this latter case, if a non negligible number of rGCs are, 
instead, misclassified bGCs, then the overall rGG average 


spectroscopic metallicity would be higher. This would hap¬ 
pen if the rGG and bGC colour distributions are not sym¬ 
metrical, with the distribution of bGCs showing a significant 
extension to redder colours. However, we note that this ex¬ 
planation is purely speculative, since the only models show- 
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6 Pastorello et al. 


ing a possible asymmetry between the colour distributions of 
the two subpopulations we are aware of show a possible con¬ 


tamination of bGCs by rGCs and not the opposite (Forte, 


Faifer & Geisler 20071. In any case, the offset between photo¬ 


metric and spectroscopic rGG metallicities is comparable to 
the uncertainties on the actual spectroscopic measurements 
an the systematic uncertainties in the colour-metallicity con¬ 
versions. We conclude that, at least in the galaxies for which 
we have a high number of spectroscopically measured GC 
metallicities, the spectral subsample is a reasonable repre¬ 
sentation of the true metallicity distribution. 


3 RESULTS 

In order to disentangle the different behaviours of the stellar 
and GC metallicities with galaxy mass, we split our sample 
into two different stellar mass bins. The radial metallicity 
profiles for the stellar and GC components of the 12 galaxies 
in this work are presented in Figure]^ (for the galaxies with 
M* > 10^^ Mq) and Figure I (for the galaxies with M* 
10 “ Me). 

Our [^/H] stellar radial proHles are extracted from the 
kriging maps presented in Pastorello et al. (20141. In ad¬ 


dition, for 7 galaxies we have obtained their inner radial 
metallicity profiles from SAURON/ATLAS®® maps. While 
the former do not extend to the galaxy centre, the latter 
do not extend to 1 Re and are, thus, complementary to 
our data. The stellar metallicity profiles, their correction 
to match with the inner SAURON/ATLAS®® proHles and 
the stellar gradient trends with galaxy mass (i.e. steeper in 
lower mass galaxies) are discussed in|Pastorello et ah (20141. 


3.1 Halo stellar metallicity from resolved RGB 
stars 

Here we compare our integrated stellar metallicities and 
GC metallicities with those from resolved RGB stars from 
HST observations. A caveat for these studies is that the 
RGB measurements are usually incomplete at high metal¬ 
licities, and therefore their values may underestimate the 
mean metallicities of the real stellar population. For two of 
our galaxies, NGC 3115 and NGC 3377, RGB metallicities 
are available and we include these values in Figure]^ 


Elson (1997 \ presented the RGB stellar metallicity for 


NGC 3115 measured from an HST/'WFPC2 field at i? « 
25 Ro. She found two distinct stellar populations, with iron 
abundances of [Fe/H] ~ —0.7 and —1.3 dex, respectively. In 
order to plot her results in our Figure [3| we convert them 
with the [Thomas, Maraston Sz Bender ( |2003 1 relation: 

l^/H] = [Fe/H] -b 0.93 [a/Fe] (3) 


where [a/Fe] is the a-element abundance. Norris, Sharpies] 
fc Kuntschnerj ( [2006 1 measured a radially constant [a/Fe] « 
0.17 dex along the NGC 3115 major axis, while they found 
a variable a-element abundance along the minor axis in the 
range 0.17 ^ [a/Fe] ^ 0.30 dex. To deal with this variabil¬ 
ity, we assume [a/Fe] = 0.24 dex. We note that the low- 


metallicity subpopulation presented by Elson (19971 has a 


comparable number of stars to the high-metallicity subpop¬ 
ulation. This may be due to instrumental calibration issues. 


More recently, for the same galaxy. Peacock et al.[ ( [2015[ ) 
obtained RGB metallicities for three different HST/ACS 
and HST /WFPC3 fields. These are also plotted in Figure]^ 
We confirm that the bulk of their measurements are consis¬ 
tent with a metallicity in between the red and the blue GC 
metallicities at the same radii. In addition, the low metal¬ 
licity tails in the three fields extend into the bGC metal¬ 


licity domain. Peacock et al. (20151 also found evidence for 


a low-metallicity stellar population in their outermost peak 
(at R = 53.8 Re), with an average metallicity value (i.e. 


yZ/H) « —1.2 dex) similar to that of Elson (19971, albeit 
with the number of RGB stars significantly lower than that 
presented in Elson ( 1997[ ). Interestingly, this low-metallicity 
stellar population has an average metallicity very similar to 
that of the bGCs (i.e. yZ/H) = —1.27 dex). 

Similarly, we add [Harris et al. ( 2007a| RGB stellar pho¬ 
tometric metallicities from HST/ACS imaging in the outer 
regions of NGC 3377. The targeted three contiguous fields, 
all yielding to a consistent metallicity of yZ/H\ « —0.6 dex. 
This is in conflict with the general trend of our SKiMS 
metallicities up to 2 Rc, which instead show a steep decline 


down to yZ/\T\ ~ —1.0 dex. Since the Harris et al. (2007al 


helds probed only a small fraction of the halo, a possible az¬ 
imuthal variation or substructure of the RGB stellar metal¬ 
licity may explain the observed difference. 

3.2 GC spectroscopic metallicity gradients 

In Figure]^ we show the combined metallicity profiles for 
both the stellar and the GC components for our two galaxy 
stellar mass bins. To obtain the combined stellar metallicity 
profiles in the two mass bins, we stack the stellar metallicity 
profiles of the single galaxies and derive the average pro- 
Hle (and the standard deviation) within logarithmic radial 
bins of 0.05 dex. The combined rGG (bGC) metallicity pro¬ 
Hles are obtained stacking the metallicities of all the rGCs 
(bGCs) in each galaxy mass bin. We derive the average pro¬ 
file (and the standard deviation) within radial bins with a 
constant number of elements (i.e. Abin = 20). 

In Tablej^we present the metallicity gradients measured 
for the rGG and bGC subpopulations in the two galaxy 
mass bins with separation at M* = 10®® Me- In the low 
mass galaxy bin, the rGG and bGC gradients are consistent 
within the associated uncertainties. Likewise, the rGC and 
bGC metallicity gradients in the high mass bin are consistent 
with each other within the associated uncertainties, and are 
shallower than the rGC and bGC metallicity gradients in the 
low mass bin. In particular, both rGC and bGC gradients 
show a strong difference between the two mass bins (i.e. 1.7a 
and 2.3 ct, respectively). Such gradients are measured in the 
whole available radial range for the different GC subpopula¬ 
tions (we find that using a restricted, common radial range 
does not change the conclusions). 

It is useful to remember that, in general, 3D metallicity 
profiles are ~ 10% steeper than those measured for projected 
galactocentric radii (Liu et al.|2011 1. 


3.3 Comparison with metallicity gradients from 
simulations 

From hydrodynamical high-resolution, cosmological zoom 


as later pointed out by Kundu & Whitmore (19981. 


simulations, Hirschmann et al. (20151 obtained stellar metal 
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Figure 2. Metallicity profiles for the high-mass galaxies in our sample. The panels (A-F) show metallicity versus elliptical galactocentric 
radius scaled by Re. The SKIMS stellar metallicity profiles are presented as green solid lines. When available, the SAURON stellar 
metallicity radial profiles are presented as a magenta dashed line. For NGC 4365 (panel D), the magenta dashed line shows instead the 
ATLAS®^ stellar metallicity profile. Blue crosses and red circles show, respectively, blue and red GC metallicities. The blue crosses and 
red circle points with errorbars show the average uncertainties on the metallicities of blue and red GCs, respectively. The red and the 
blue dotted lines show the average metallicities of the rGC and bGC subpopulations, respectively. 



Low-mass V[Z/H] 

High-mass V[Z/H] 


(dex/dex) 

(dex/dex) 

rGC 

-0.26 ±0.08 

-0.07 ± 0.08 

bGC 

-0.45 ±0.13 

-0.08 ± 0.09 


Table 2. GC metallicity gradients in the two galaxy mass bins 
with separation at M* = 10^^ Mq. The gradients are measured 
in the whole available radial range (i.e. 0.6 R 14 Ro) 


licity profiles out to large radii for 10 galaxies with stellar 
mass in the range 1.7 x lO^'^ < M* < 3.7 x 10^^ M©. They 
found that when including stellar wind feedback the average 
stellar metallicity profile slope is consistent with the photo¬ 
metric observations by La Barbera et al. (2012|. Hence, we 


compare our metallicity gradient results with the predictions 
from only their simulations including stellar wind feedback. 
In what follows we assume that the |Hirschmann et al.| ( [^15[ ) 
gradient predictions are predominately driven by the differ¬ 
ent accretion histories of the simulated galaxies, and not by 
the mass (which has a very narrow range). 


We do not compare the absolute metallicity values of 


Hirschmann et al. (20151, although we do apply an offset in 


order to have the profiles match photometric observations 
from the literature. This is because the absolute value of 
stellar metallicity in simulations is sensitive to the subgrid 
baryonic physics. Therefore, it is not safe to directly compare 
our absolute values with those from simulations, while in- 
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log (i?/Rc) [dex] log (i?/Rc) [dex] 


Figure 3. Metallicity profiles as in Figure]^ but for the 6 low-mass galaxies in the sample. For NGC 3115 (panel B) and NGC 3377 
(panel C) the density plots show the photometric stellar metallicity distribution functions (MDF) from [Peacock et al.| ( |201^ and |Harris| 
|et al.| | [2007a[ |, respectively. The radial position of the density plots corresponds to the centre of the observed fields. For better readability, 
in both cases the radial extension of the density plots does not represent the actual size of the observed field. These density plots are 
colour-coded according to the colour bars in the same panels. For NGC 3115 (panel B) we also show the photometric metallicities of the 
halo metal-rich and metal-poor stellar populations as found by |Elson| l|1997^ as red and blue triangles, respectively. 


stead we can qualitatively compare the metallicity gradients 
(which depend mostly on the galaxy formation history). 

( 2015|), the average 


history from the Hirschmann et al. (20151 simulations with 


In figure 4 of Hirschmann et al. 


metallicity profile is arbitrarily scaled in metallicity for vi¬ 
sualization purposes (M. Hirschmann, private communica¬ 
tion). We extract from this figure the average metallicity 
profile and offset it to match the observed prohle for the 
high-mass galaxy sample in figur e 8 of|La Barbera et al.| 
(20121. In particular, we adopt the La Barbera et al. (20121 
metallicity profile obtained with the [Bruzual fc Chariot] 
(20031 single stellar population models. In this way we ob¬ 


tain the vertical offset to apply to the [Hirschmann et ah] 
| |2015| ) prohles. 

In Figure]^ we plot two extreme examples of assembly 


stellar wind feedback at z = 0, representative of minor and 
major merger dominated evolutionary history, respectively. 
In particular, in the left panel (i.e. our low-mass galaxy bin) 
we plot the stellar metallicity profile from the minor merger 
model M215, which is the model with the highest fraction 
of minor mergers and no past major merger event. In the 
right panel (i.e. our high-mass galaxy bin) we plot the stellar 
metallicity profile from the major merger model M305. This 
is the model with the highest accreted mass fraction from 
major mergers. In particular, it includes a major merger 
event at ^ « 0. In both cases, the vertical offsets applied 
to the simulated stellar metallicity profiles provide a good 
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M* < 10“ Mq M* > 10“ Mq 



il/R, R/R, 


Figure 4. Combined stellar, rGC and bGC radial metallicity profiles for the low-mass {left panel) and high-mass {right panel) galaxies 
in our sample. In both panels, the combined stellar, rGC and bGC metallicity profiles are presented as a green dotted, a red dashed 
and a blue solid line, respectively. The green curves have been obtained combining SKiMS and SAURON/ATLAS^^ metallicity profiles. 
The green, red and blue shaded region represent the la confidence limits in the stellar, rGC and bGC cumulated metallicity profiles, 
respectively. In the stellar case, the confidence intervals are dominated by the scatter among the galaxies. In the rGC and bGC cases, the 
confidence intervals are dominated by the scatter among the GC metallicities within the same galaxy. The black dotted line in the left 
panel and the black dashed line in the right panel show the stellar metallicity profiles for the model M215 and M305 of |Hirschmann et al.| 
l |2015[ |, respectively. Such models represent two extreme formation histories (i.e. M215 experienced only minor mergers while M305 has 
had a recent major merger). These simulation stellar metallicity profiles have been scaled in [Z/H] to match with the photometrically 
observed metallicity profiles of the high-mass sample in |La Barbera et al. j poT^ . 


match with our spectroscopically measured stellar metallic¬ 
ity profiles at the common radius. 


Hirschmann et al. (20151 provided the stellar metal¬ 


licity gradients in the range 2 < i? < 6 Re. The stellar 
metallicity profile of M215 has a gradient of —0.41 dex/dex, 
consistent with that of bGCs in our low-mass bin (Table [^. 
The stellar metallicity gradient of M305 has a gradient of 
—0.13 dex/dex, consistent with both the gradients of rGCs 
and bGGs in our high-mass bin (Table |^. 

The range of metallicity gradients in the minor merger 
simulations by [Hirschmann et al.| | |2015| ) is —0.22 ^ 

V[.Z/H] ^ —0.76 dex/dex, while the range in the major 
merger simulations is —0.10 V[.^/H] ^ —0.46 dex/dex. 

Within these ranges, both our high-mass bin rGG and bGC 
metallicity gradients are fully consistent with those from the 
major merger simulations. However, in the low-mass bin our 
GC subpopulation metallicity gradients are consistent with 
both major and minor merger simulations. We will return 
to these comparisons in more detail in Section [4.3| 


3.4 The metallicity offset between inner stars and 
rGCs 

Thanks to our radially extended stellar and GC spectro¬ 
scopic metallicity sample we are able to measure the metal¬ 
licity offset at 1 Re between stars and rGCs. 

We verify the existence of a correlation between the 
stellar-rGC metallicity difference with galaxy mass from our 
sample of 12 galaxies. In Figurej^we plot the metallicity dif¬ 
ference A [Z/H] between the rGCs and the stars at 1 Rc for 


Galaxy 

m 

<? 

(dex) 

NGC 1407 

-0.01 ± 0.02 

-0.31 ± 0.08 

NGC 3115 

0.00 ±0.01 

-0.18 ±0.10 

NGC 4278 

-0.02 ±0.01 

-0.40 ±0.08 

NGC 4365 

-0.02 ±0.02 

-0.37 ± 0.07 

NGC 4649 

0.03 ±0.04 

-0.38 ±0.12 


Table 3. Best linear fit parameters for the rGC metallicity pro¬ 
files of the five galaxies with the highest number of rGCs. The 
linear fit is in the form: [Z/H] = mlog(R/Re) -|- q 


our sample as a function of different galaxy mass proxies. 
In particular, for the five galaxies with a high number of 
rGCs (i.e. NGC 1407, NGC 3115, NGC 4278, NGC 4365 
and NGC 4649), we have been able to extrapolate the rGC 
metallicity to R ~ 1 Rc from the best linear fit of the 
rGC metallicity profile. For such galaxies, the best linear 
fit parameters are presented in Table Because of the low- 
number of rGCs in the other galaxies, we adopt the average 
metallicity of the whole rGC subpopulation in each galaxy 
as a proxy for the value at R ~ 1 Ro. 

In the first panel we use galaxy central velocity disper¬ 
sion (To from Table The best fitting linear function is: 

A[Z/H] = (2.50 ± 0.75) log(To - (5.71 ± 1.76). (4) 

The rms of this relation is 0.08 and the Spearman index 
value (a statistical measure of the correlation between two 
variables) is rg = 0.45 with a significance of p = 86.2 per 
cent. 
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Figure 5. Difference between stellar and rGC metallicities (i.e. [Z/i?]star — [-^/^IrGc) versus various galaxy mass proxies. The stellar 
metallicities values are obtained at 1 Re. For 5 galaxies we have been able to extrapolate the rGC metallicity to 1 Re from a linear 
fit of their radial metallicity distribution. These galaxies are presented as red points. For the other galaxies, the rGC metallicities are 
obtained as the average metallicity of the whole rGC system. These galaxies are presented as black points. In the left panel the metallicity 
differences are plotted against the stellar central velocity dispersion, in the central panel against the total stellar mass and in the right 
panel against the total B-band magnitude of the host galaxy. The dashed lines are the best-fitting linear functions in the three cases. A 
mild correlation between the stellar-rGC metallicity difference and galaxy mass is seen. 


In the second panel we use total stellar mass M* from 
Table obtaining the best-fit linear relation: 

A[Z/H] = (0.50 ± 0.32) log M* - (5.42 ± 3.50). (5) 

The rms is 0.24 and the Spearman index is rs = 0.41 with 
a significance of p = 80.9 per cent. 

Lastly, in the third panel we use total B-band magni¬ 
tude of the host galaxy, as presented in Table The best-fit 
relation is: 


A[Z/H] = (-0.19 ± 0.14) logMn - (3.84 ± 2.81). (6) 


The rms is 0.28 and the Spearman index is rs = —0.37 with 
a significance of p = 77.0 per cent. 

In summary, we find that the difference between the 
stellar metallicity and the rGC subpopulation metallicity 
at 1 Re correlates mildly with the galaxy mass. Thus, the 
rGC subpopulations in low-mass galaxies show an average 
metallicity closer to that of the galaxy stars at similar radii 
than in high-mass galaxies. 

This trend may suggest a longer lag between the peaks 
of formation for rGCs and stars in more massive galaxies 
(i.e. continued enrichment in galaxy stars), as suggested by 


Spitler (20101. Supporting this view, Montes et al. (20141 


found that the rGCs in M87 are older and more metal-poor 
than the bulk of the stars in the central regions of the galaxy. 
However, as lower-mass galaxies have been found to expe¬ 
rience a more extended star formation than their higher- 


mass counterparts (Thomas et al. 2005 McDermid et al. 


2015), the trend between stellar-rGC metallicity difference 


and galaxy mass should be the opposite (i.e. larger metal¬ 
licity offsets in lower mass galaxies). 


However, there is another possible cause for this corre¬ 
lation. The correction applied to our stellar metallicities is 
obtained from the offsets with SAURON metallicities (see 


Section [2^. Since this offset is probably linked with a more 
bottom-heavy IMF in the central regions of high-mass galax¬ 
ies (Pastorello et al. 2014), this may not apply to the re¬ 
gion at R ^ 1 Re, where the IMF has been found to be 
shallower and similar among galaxies with different mass 
I Martin-Navarro et al.|2015 |. In this case, the stellar metal¬ 
licity at 1 Re in high-mass galaxies may be systematically 
overestimated. 


Goudfrooij & Kruijssen (2013) found that a colour offset 


between rGCs and stars in high-mass galaxies exists out 
to around 1 Ro. From the Lick indices analysis in two of 
their galaxies, they found consistent ages and metallicities 
between rGCs and stars. Therefore, they claimed that only 
a steep IMF for the galaxy stars can explain the colour offset 
with the rGCs. This is in conflict with our findings at 1 Re, 
since we spectroscopically measure different metallicities for 
rGCs and stars. However, as discussed above, a shallower 
IMF slope in high-mass galaxies at 1 Re may reduce the 
stellar metallicity at such radii. 


3.5 Summary of the results 

In summary, we find that the possible metal-poor stellar 
population at large radii from RGB stellar studies shows 
a metallicity similar to that of bGCs. In addition, we ob¬ 
serve that the metallicity gradients of bGCs and rGCs are 
consistent one with another within the uncertainties in the 
same galaxy mass bin. Moreover, in the high-mass galaxy 
bin, both bGCs and rGCs show shallower metallicity gradi¬ 
ents than in the low-mass galaxy bin. We also compare our 
metallicity profiles in low- and high-mass galaxies with stel¬ 
lar metallicity profiles from simulations assuming a minor 
merger and a major merger dominated evolutions. We find 
that the bGC metallicity gradient in our low-mass galaxies 
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is consistent with the former, while in the high-mass galax¬ 
ies both rGC and bGC metallicity gradients are consistent 
with the latter. Moreover, we observe that the average stellar 
metallicity at 1 Re is very similar to that of the rGC sub¬ 
population, with larger differences for higher-mass galaxies, 
where the metallicity of the stars is higher than that of the 
rGCs. In the few cases where we can probe the stellar com¬ 
ponent to large Ro (i.e. NGC 3115, NGC 3377, NGC 4473 
and NGC 4494), we find that the stellar metallicity at large 
radii tends to be similar to the average metallicity of the 
bGC subpopulation. 


4 DISCUSSION 

4.1 Stellar metallicity profiles at large radii 

In what follows, we discuss the results within the frame¬ 
work of the two-phase galaxy formation scenario, although 
the general conclusions we draw from our observations can 
be used as constraints for future different galaxy formation 
models. 

In a two-phase formation scenario, the outer regions of 
massive ETGs are thought to be dominated by stars formed 
ex-situ, i.e. accreted from smaller companion galaxies or dur¬ 
ing major merger events. In the first case, the relation be¬ 
tween a galaxy’s mass and its mean stellar metallicity dic¬ 
tates that small accreted galaxies are made up of metal-poor 
stars. Therefore, if the stars in the outskirts of ETGs are 
mostly accreted from smaller satellite galaxies, these regions 
should reveal a stellar population with a low mean metallic¬ 
ity. In this case, higher metallicity inner regions and lower 
metallicity outskirts in the final galaxy are expected to be 
separated by a transition region. The metallicity profile in 
this region will be steeper than the original one at the same 


radii (Hirschmann et al. 20151. However, quantitative pre¬ 


dictions about the galactocentric radius of this region vary 
in the literature. One issue is the estimation of the effects of 
SNe and AGN feedback on the final stellar metallicity pro¬ 
files, even though some recent studies have started to tackle 
this problem in detail (e.g. [Hirschmann et al.|2015p . 

In the case of major mergers, accreted stars will have 
a metallicity comparable with that of pre-existing in-situ 
formed stars. Thus, any pre-existing metallicity gradient will 
be washed out by the major merger event ( [Hirschmann et ah] 
2015). In addition, a significant fraction of in-situ stars is 


pushed at large radii, further diluting their metallicity gra¬ 
dient. 

In [Pastorello et alT] ( [2014[ ) we found that the stellar 
metallicity outside 1 Rc declines steeply with radius in low- 
mass galaxies, while it displays a shallower gradient in high- 
mass systems. This is either a reflection of a relative lack 
of mergers in the assembly history of our lower-mass ETGs 
and/or we are probing the region between in-situ and ex-situ 
dominated regions, where the latter are built-up via minor 
mergers. If it is a lack of mergers, then the galaxy metallicity 
profiles are still similar to those formed during the dissipa¬ 
tive phase of galaxy formation, a If, instead, minor mergers 
dominate in the evolutionary history of low-mass galaxies, 
then a stellar low-metallicity “plateau” should be present at 
larger radii, where the accreted stars dominate the stellar 
population. 


Spectroscopically probing such regions is still a chal¬ 
lenge, but some hints of the stellar metallicity profiles at 
large radii come from resolved RGB stars in nearby low- 
mass galaxies, albeit with some limitations (e.g. limited ra¬ 
dius baseline, observational completeness issues). For exam¬ 
ple, the almost radially constant RGB stars’ metallicities in 


NGC 3115 (Elson 1997 Peacock et al. 2015) suggest that 


the average integrated stellar metallicity in these outer re¬ 
gions does indeed reach a “plateau” for R > IR^. The RGB 
metallicities obtained for NGC 3377 at large radii by [Har-j 
ris et al. (2007a I also show a flat metallicity radial profile. 
On the other hand, in the nearby galaxy NGC 5128, RGB 
metallicities have been found to show a mild gradient (i.e. 
A[.Z/H]/ARe « —0.030±0.003 dex) out to more than 20 Ro 
( [Rejkuba et al.|2014 Bird et 5/][2015 ). 

An interesting addition to this picture is the observa¬ 
tion of distinct stellar components in the data of both jElsonj 
(1997) and Peacock et al. (2015). Their median metallicities 


are qualitatively close to those of the red and blue GCs in 
NGC 3115 at the same radii. 


4.2 What is the origin of GC metallicity 
bimodality? 

Most GC systems of large galaxies show a bimodality in 
colour, which generally reflects a bimodality in metallicity 


I Brodie fc Strader|2006 Brodie et al.|2012 Usher et al.|2012' 
and references therein). Most current theoretical models of 
GC formation predict that rGCs are mostly formed in-situ, 
together with the bulk of the host galaxy’s stars, and bGCs 
form in dwarf galaxies that are later accreted, preferentially 


in the outer regions of the larger galaxies (Cote, Marzke 


fc West||1998 Tonini|[2013 Katz fc Ricotti|[2013 2014), al 

though some in-situ GC formation in the galaxy halo may 
occur ( Canning et al.|2014 |. In this scenario, the bGC aver¬ 
age metallicity is lower than that of rGCs because it reflects 
that of the low-mass systems where the bGCs form. 

However, another possible scenario for the bimodality 
formation allows both rGCs and bGCs to be formed in-situ 
and also accreted from low-mass galaxies. In this case, in- 
situ rGCs and bGCs mostly populate the inner regions of 
the host galaxy, while the accreted GCs lie in the outer re¬ 
gions. If the main sources of ex-situ GCs are minor mergers, 
one may expect a lower number of accreted rGCs than ac¬ 
creted bGCs. This because low-mass accreted galaxies likely 
had low numbers of rGCs, hence their contribution to the 


host galaxy’s rGC subpopulation should be small (Forbes, 


Brodie fc Grillmalr||1997| [Strader et al.|[2006| [Peng et al. 
2006[ Fbrbes et al.pOll ). On the other hand, major merg¬ 
ers will contribute with a more similar number of rGCs and 
bGCs. 

The most notable difference between the two scenarios 
concerns the origin of the bGCs, for which the available the¬ 
oretical models predict almost exclusively ex-situ formation, 
while in the second case a significant fraction of bGCs may 
have formed in-situ. 

Differences in rGCs and bGCs azimuthal, kinematic and 
metallicity distributions are expected in the two scenarios. 

If bGCs are exclusively formed ex-situ and later ac¬ 
creted, their azimuthal and kinematic distributions should 
be, to first order, decoupled from those of the rGCs and the 
galaxy inner stars. Instead, if bGCs and rGCs formed under 
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similar conditions, both will show in the inner regions kine¬ 
matics and azimuthal distributions comparable to those of 
the galaxy stars (i.e. with similar ellipticity and PA). Ob¬ 
servations of GC kinematics and azimuthal distributions do 
not give a clear answer to the origin of bGCs. In fact, several 
observations favour the latter scenario (i.e. bGCs and rGCs 
formed under similar conditions), with galaxies displaying 
azimuthal distributions for both rGCs and bGCs similar to 


those of the host galaxy’s stellar component (e.g. [Dirsch, 


Schuberth & Richtler 2005 Hargis & Rhode 2014^, while 


others have more decoupled components (e.g. Strader et al. 
2011 Kartha et al.|20l4l. Moreover, in a number of nearby 


ETGs the kinematics of bGCs have been found to be signif¬ 
icantly distinct from those of the stellar component and the 
rGCs (e.g. Schuberth et al.|2010 Pota et S(]|2013 1. 

A better way to disentangle the contribution of in-situ 
and ex-situ formed GCs in the two subpopulations is to 
examine the relation between their average metallicity and 
host galaxy mass. If a GC subpopulation exclusively formed 
in-situ, then its metallicity should show a positive correla¬ 
tion with host galaxy mass, i.e. in high-mass galaxies the in- 
situ formed GCs will have higher average metallicities than 
those in low-mass galaxies. On the other hand, the average 
metallicity of an accreted component would show a weaker 
correlation with the current host galaxy mass. This is be¬ 
cause the mass of an accreted satellite is not directly depen¬ 
dent on the mass of the final host galaxy, although high mass 


licity profile in the region 2 Si i? ^ 6 Re is steeper in the 
minor merger simulations, since the accreted stars at large 
radii will have a lower metallicity than the in-situ formed 
stars at the same radii. In general, one may expect a similar 
trend for the case of in-situ and ex-situ GCs, with steeper 
gradients if these latter are accreted during minor merger 
episodes. However, the specific frequency of rGCs is lower 
in low-mass galaxies and, therefore, minor merger events 
will generally contribute with more bGCs than rGCs. Ma¬ 
jor mergers, instead, will strongly reduce the steepness of 
the metallicity profiles at almost any radius. This is because 
the accreted GCs will be radially mixed in the host galaxy 
and the in-situ formed GCs will be pushed at larger radii. 
However, since the GCs formed in more massive satellites 
tend to be accreted at smaller galactocentric radii, a mild 
metallicity gradient may be created by this mass segregation 
effect. 

Therefore, if rGCs are mostly formed in-situ and bGCs 
are formed ex-situ, the former should show steeper metal¬ 
licity gradients than the latter. Instead, if the two GC sub¬ 
populations have had a similar formation process, both be¬ 
ing formed in-situ and accreted, they are expected to show 
similar metallicity gradients over all radii (i.e. steep in the 
inner regions and shallow in the outer regions). 

In the literature there are several massive ETGs for 
which both the rGCs and bGCs in the inner regions show 
a significant colour/metallicity gradient, while in the outer 


galaxies tend to accrete higher mass satellites (Oser et al. regions such gradients flatten (Forte et al. 


2012). So, if rGCs are mostly formed in-situ and bGCs are 


accreted, the rGCs will show a strong trend of their average 
metallicity with the host galaxy mass, while for the bGCs 
this trend will be weak/negligible ( |Bekki et al. [200^ 1. 

It is not clear whether or not the bGC average metal¬ 
licity trend with host galaxy mass is consistent with 
this picture. A weak correlation between bGG average 
colour/metallicity and galaxy mass has been confirmed by 
several studies (e.g. Larsen et al. 2001[ Strader, Brodie & 
Forbes 2004| Strader et al. [2006 


2001 


Harris 


2009a|b| Forbes et al.|2011 

Arnold et al. 

201ll 

Liu et al. 

2011 

[Jennings et al.|2014||I 

Largis & Rhode 

2^41 p 

Mta et al. 


20151. While this strongly supports a common in-situ forma- 


Peng et al. 2006). How¬ 


ever, in the cases where a colour/metallicity gradient exists 


(e.g. Peng et al. 20061 this result may be simply a conse¬ 


quence of an aperture bias (i.e. a GG system average will 
lead to measuring the colours/metallicities of mostly the in¬ 
nermost bGGs in massive galaxies). Once corrected by this 
bias, the estimated slope for the relation between bGC av¬ 
erage metallicity and host galaxy mass may be greatly re¬ 
duced, potentially disappearing, as discussed in |Liu et ah] 
(2011). This would reinforce the theoretical view of exclu¬ 


sively externally formed bGCs. In our low-mass galaxy bin 
we note that bGCs show a lower average metallicity (i.e. 
\Z/}¥\ = —1.21 ±0.04 dex) than in the high-mass galaxy 
bin (i.e. [Z/Yi\ = —1.12 ± 0.03 dex) if compared over the 
same radial range (i.e. 1 ^ R/Ro ^ 8). Since our spectro¬ 
scopic sample is not affected by aperture bias, a mild trend 
of bGC average metallicity with galaxy mass may exist. 

Another key difference between in-situ and ex-situ 
formed GCs concerns their radial colour/metallicity gradi¬ 
ents. Similar to the predictions for the stellar component, in- 
situ formed GCs will be more metal-rich closer to the centre 
of the gravitational well. This is because the gas from which 
GCs form is more metal-rich where the stellar evolution is 
faster. 


tion for the two GC subpopulations in the inner regions, the 
metallicities have been obtained from photometric colours. 
Even though most GCs have old ages, these colour gradients 
may include age variations, thus obscuring the true metal¬ 
licity gradients. In summary, from the available literature 
it is not clear whether bGCs may be also formed in-situ or 
are exclusively accreted (as suggested by most theoretical 
models). 

In Section [3.2| we presented the spectroscopic metallic¬ 
ity gradients for both rGCs and bGCs in two galaxy mass 
bins. In both mass bins, the gradients of the two GC sub¬ 
populations are similar to each other. These results support 
a scenario in which rGCs and bGCs in the same galaxy 
share a common origin. In particular, the significant bGC 
and rGC metallicity gradients in our low-mass galaxy bin 
are in agreement with a mostly in-situ origin for both the 
GC subpopulations. 

4.3 Fraction of accreted GCs 

As seen in the previous Section, from the rGC and bGC 
metallicity gradients it is possible to infer the existence of 
both in-situ and ex-situ formation for both GC subpopula¬ 
tions. Under the two-phase galaxy formation scenario, the 
ex-situ fraction of stars is predicted to be higher in more 
massive galaxies than in their low-mass counterparts (e.g. 
Oser et al.||2010 Hirschmann et al.|2013 |. The Hirschmann 
et al. (2015) simulations, including stellar wind feedback. 


Hirschmann et al. (2015) found that the stellar metal- 


also showed that the ex-situ stars accreted during minor 
merger episodes contribute to steepening the pre-existing 
in-situ stellar metallicity gradient. In addition, major merg- 
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ers create a shallower stellar metallicity profile in the outer 
regions. If these trends with the formation history apply 
also to the GC subpopulations, one can draw some simplis¬ 
tic conclusions from the GC metallicity gradients in galax¬ 
ies of different mass. In particular, if we assume that GCs 
are a good proxy for the stellar component at large radii, 


of the minor and major merger models. We note that, even 


then we can compare the Hirschmann et al. (20151 stellar 


metallicity gradient predictions with our GC metallicity gra¬ 
dient. A caveat to keep in mind is that [Hirschmann et ^ 


(20151 simulations regard galaxies with M > IO'^^Mq, and 


thus our low-mass galaxies are not directly comparable with 
their simulations. However, in what follows we assume that 


Hirschmann et al. (20151 predictions are also valid for the 


galaxies in our low-mass bin (see Section 3.31. 

Similarly to the stellar component, if in high-mass 
galaxies a high fraction of rGCs and bGCs are accreted 
during major mergers, then the final rGC and bCC metal¬ 
licity profiles should be fiat. Conversely, if most of the GCs 
form in-situ or are accreted during minor mergers, then they 
should show steep metallicity gradients. In particular, if the 
ex-situ GCs are accreted during minor mergers episodes, 
then the bGC metallicity profiles may be expected to be 
steeper than those of the rGCs in the same galaxies. This is 
because, given the different specific frequency of rGCs and 
bGCs in low-mass systems, minor merger episodes will in 
general contribute with more accreted bGCs than rGCs. As 
a consequence, the total metallicity gradients of both GC 
subpopulations will be steeper than in the major merger ac¬ 
cretion case, but rGCs will also show a shallower gradient 
than the bGCs. 

This is consistent with our results from Section r3.2l The 
shallower (consistent with being flat) metallicity gradients 
we measure for the rGCs and bGCs in high-mass galaxies are 
consistent with the presence of a higher fraction of ex-situ 
formed GCs accreted during major mergers in these systems. 
We find that both these metallicity gradients are consistent 
with the stellar metallicity gradients obtained in a major 
merger dominated evolutionary history in the simulation by 


Hirschmann et al. (20151. 


On the other hand, the steeper metallicity gradients 
for rGCs and bGCs in low-mass galaxies show consistency 
with both minor and major merger dominated formation 
histories. This is qualitatively similar to that observed for 


the stellar component in low-mass galaxies (Pastorello et al. 
|2014| ). The rGC metallicity gradient in the same galaxy mass 
bin is shallower. We argue that this is because the number 
of accreted rGCs is lower than for the bGCs and, therefore, 
their metallicity gradient is more similar to the in-situ stellar 
metallicity gradient. 

In Figure we present a further comparison of our GC 
metallicity gradients with the simulations by jHirschmannj 
et al.|([2015[) that include stellar wind feedback. In particular, 


we plot the simulated gradients for both of the minor and 
major merger models against the mass weighted merger ratio 
and the fraction of ex-situ stars in the simulations. 

In general, the simulated minor merger stellar metallic¬ 
ity gradients have lower mass weighted merger ratios and 
lower fractions of ex-situ stars than the simulated major 
merger models. In the same Figure, we also show the aver¬ 
age metallicity gradients of rGCs and bGCs in our two mass 
bins. In the low-mass galaxy bin, the average metallicity for 
the two GC subpopulations is qualitatively in between those 


though the Hirschmann et al. (20151 simulations did not in¬ 


clude galaxies in our low mass bin mass range, the effects 
of the formation history on the metallicity gradients is the 
dominant factor. In the high mass galaxy bin, the average 
metallicity gradients of both the two GC subpopulations are 
similar to those from the major merger models. 


5 CONCLUSIONS 

Here we present for the first time a systematic compilation of 
spectroscopically measured metallicities for the stellar and 
the globular cluster components of 12 galaxies as part of the 
SLUGGS survey. Metallicities at large radii for both stars 
and globular clusters have been obtained from the CaT in¬ 
dex. For the in nermost radii, we a dopt stellar met allicities 
from SAURON ( Bacon et al.j2001 | and ATLAS^° (Cappel- 
lari et al. surveys for the galaxies in common. Com¬ 

bining SAURON, ATLAS®° and SLUGGS values, we have 
been able to spectroscopically measure the metallicity of the 
stellar component out to almost 4 Ro, and of the red and 
blue globular clusters out to almost 15 Rc. 

Our main results can be summarized as follows: 


• We spectroscopically confirmed that some red globular 
clusters exist out to i? > 10 Rc. 

• We compare the average stellar metallicity at 1 Rc with 
that of the red globular cluster component, finding that 
high-mass galaxies host a stellar population that is more 
metal-enriched than the red globular cluster component, 
while in low-mass galaxies the stellar and the red globu¬ 
lar cluster metallicities are more similar. We claim that this 
can not be simply due to a bottom-heavy IMF. 

• We find that the metallicity gradients of red and blue 
globular clusters are consistent within the same galaxy mass 
bin. In particular, in the low galaxy mass bin blue and red 
globular cluster metallicity gradients are —0.26 ± 0.08 and 
—0.45 ±0.13 dex/dex, respectively. In the high galaxy mass 
bin, red and blue globular cluster metallicity gradients are 
—0.07 ± 0.08 and —0.08 ± 0.09 dex/dex, respectively. This 
consistency between globular cluster subpopulation metal¬ 
licity gradients within the same galaxy mass bin suggests 
a common formation process for both red and blue glob¬ 
ular clusters in the same systems. The colour/metallicity 
bimodality is then potentially a consequence of different 
epochs for the two subpopulations, as already concluded by 


photometric studies (e.g. Arnold et al. 2011 Forbes et al. 


20111. This is in contrast with most of the current theo¬ 


retical models that consider red globular clusters as mostly 
formed in-situ and blue globular clusters as mostly formed 
ex-situ and later accreted into their present day host galaxy. 

• We measure metallicity gradients for both the GC sub¬ 
populations that are steeper in the low galaxy mass bin than 
in the high galaxy mass bin. From the comparison with stel¬ 
lar metallicity gradients from simulations at similar radii, we 
find that this trend is qualitatively consistent with a mass- 
dependent galaxy formation history, in which the ex-situ 
fraction and the importance of major mergers increase with 
the galaxy mass. In particular, in high-mass galaxies the 
higher number of accreted globular clusters yields shallower 
metallicity gradients because of the dynamical mixing of in- 
situ and ex-situ formed globular clusters. On the other hand. 
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Figure 6. Outer metallicity gradients against mass weighted merger ratio and ex-situ stellar fraction. In the left panel, the stellar 
metallicity gradients for |Hirschmann et al.| ) f2015| l simulated galaxies in the wind feedback models at 2 = 0 are plotted against the mass 
weighted merger ratio. These gradients are measured in the radial range 2 ^ i? ^ 6 Re. In the right panel, the same stellar metallicity 
gradients are plotted against the total fraction of ex-situ stars in the simulations. In both panels the red circles are for minor merger 
dominated formation and blue triangles for major merger dominated formation. The solid and dashed black lines show the average 
metallicity gradients for rGCs and bGCs obtained in this work for the low and high galaxy mass bins, respectively. These gradients are 
obtained in the whole available radial range for our GCs. Qualitatively, the GC metallicity gradient in our high mass galaxies is similar 
to those in simulated galaxies with a high fraction of ex-situ stars accreted via major mergers. Our low mass galaxy GC metallicity 
gradient is similar to those in simulated galaxies with a lower fraction of ex-situ stars and more influence from minor mergers. 


the blue globular cluster metallicity gradient is consistent 
with the steep stellar metallicity gradient from simulations 
of minor merger dominated galaxy formation, while the red 
globular cluster metallicity gradient in the same galaxies is 
qualitatively shallower. This is in agreement with a galaxy 
formation dominated by minor mergers, since the merger of 
lower-mass satellites will accrete more blue than red globular 
clusters. 
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Name 

R.A. 

Dec. 

{9 - i) 

i 

GaT 

[Z/H] 


[deg] 

[deg] 

[mag] 

[mag] 

A 

[dex] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

NGC1023_GG1 

40.1543 

39.0783 
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7 

'•’^^-1.01 
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NGC1023_GG2 
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1.07 ± 0.01 

19.71 ± 0.01 

6 

—0 76+0-10 

NGC1023_GG3 

40.1460 

39.0370 

1.12 ±0.06 

22.72 ±0.04 

8 

0 O'j+O-^l 

U-UO_o 82 


Table 4. GC CaT and metallicity measurements summary. The full version of this table is provided in a machine readable form in the 
online Supporting Information. The columns present: (1) Globular cluster IDs. (2) and (3) Right ascension, declination in the J2000.0 
epoch. (4) [g — i) colour. (5) i magnitude. (6) CaT index measurement. (7) Total metallicity. 
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APPENDIX A: UPDATED STELLAR 
METALLICITY MAPS AND RADIAL PROFILES 


Since the publication of Pastorello et al. (20141, we have 
obtained new stellar metallicity data from the ongoing 
SLUGGS survey. In particular, for the purposes of this work 
we have updated the metallicity dataset for NGC 3115. In 
Figure |A1| we present the new metallicity profiles and 2D 
map. An important difference between this and the |Pas-| 


torello et al. (2014l metallicity maps for NGC 3115 is that 


the higher metallicity substructure previously extracted in 
the North-East region region of the galaxy is not visible 
anymore. The relatively low number of datapoints in the 
old version of the 2D kriging map biased the results, in the 
sense that few high metallicity datapoints in the same region 
caused an overestimation of the kriging extracted metallic¬ 
ity. 


log{R/RJ 




Figure Al. ID metallicity profiles and 2D metallicity map of 
NGC 3115. The upper panel shows the metallicity profile ex¬ 
tracted from the kriging map (black solid line) versus the galacto- 
centric radius in logarithmic space. The central panel shows the 
same profile, but on a linear radial scale. In both these panels, 
the measured data points are shown as squares colour coded ac¬ 
cording to their metallicity. The black dotted vertical line shows 
the radius corresponding to 1 Re. The blue dashed line shows 
the metallicity profile along the major axis as measured by |Nor-| 
|ris. Sharpies fc Kuntschner| l |200^ . The lower panel shows the 2D 
metallicity map from kriging and the measured data points colour 
coded according to their metallicity values. The size of each point 
is inversely proportional to its uncertainty. The black dashed lines 
show the isophotes at 1, 2 and 3Re, with ellipticity and PA from 
table 2 in [Brodie et al.] ( |2014[ |. The data points, radial profiles 
and 2D metallicity map have been corrected with the empirical 
relation presented in equation 7 ofjPastorello et al.|(|201^. 
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